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ABSTRACT 

We present a Bayesian analysis of large-scale structure (LSS) and cosmic microwave back- 
ground (CMB) data to constrain the form of the primordial power spectrum. We have extended 
the usual presumption of a scale invariant spectrum to include: (i) a running spectral index; 

(ii) a broken spectrum arising perhaps from an interruption of the potential driving inflation; 

(iii) a large scale cutoff in power as the first year WMAP results appear to indicate; (iv) a 
reconstruction of the spectrum in eight bins in wavenumber; and (v) a spectrum resulting 
from a cosmological model proposed by Lasenby & Doran, which naturally exhibits an ex- 
ponential drop in power on very large scales. The result of our complete Bayesian analysis 
includes not only the posterior probability distribution from which parameter estimates are 
inferred but also the Bayesian evidence. This evidence value is greater for a model with fewer 
parameters unless a more complicated model provides a significantly better fit to the data, 
thus allowing a powerful method of model selection. We find that those models exhibiting 
any form of cutoff in power on large scales consistently produce higher evidences than ei- 
ther the Harrison-Zel'dovich or single spectral index spectra. In particular, within the best-fit 
concordance cosmology, we find the Lasenby & Doran spectrum to show significantly larger 
evidence as compared to the other models. 

Key words: cosmological parameters - cosmology:observations - cosmology:theory - cos- 
mic microwave background - large-scale structure 
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1 INTRODUCTION 

Recent cosmological surveys of the cosmic microwave background 
(CMB), in particular observations by the Wilkinson Microwave 
Anisotropy Probe (WMAP), and large-scale structure (LSS) data, 
such as from the two degree field galaxy redshift survey (2dFGRS), 
have provided a uniqu e tool in constraining the structure an d con- 
tents of the Universe dSpergel et al.ll2003l : lRebolo et al]2004l) . One 
of the most exciting results to emerge from these analyses is pos- 
sible structure in the primordial power spectrum, features in which 
can be used to constrain early Universe physics, including infla- 
tionary models. 

The simplest model of inflationary perturbation generation in- 
volves a single, slowly rolling inflaton field producing a nearly 
scale-invariant spectrum, although even this is only true in the 
limit of an infinitely slo w roll with infinite Hubble expansion 
damping dDodelsorJ f2002). More gener ally, phase transition s from 
force unification in the early Univer se dBarriga ft alj|200lh . mul- 
tiple field inflation, hybrid inflat i on dLinddl 19831) and other mod- 
els dAdams. Ross & Sa rkar 1997; Wang 1994) all produce features 
in the primordial spectrum. It should also be noted that a scale- 
invariant spectrum does not uniquely identify inflation as the source 
of density perturbations, but a spectrum containing specific features 
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would be indicative of a specific model of inflation. Hence the study 
of the form of the spectrum is crucial in analysis of future data. 

Generalisations of the shape of the primordial spectrum 
have been attempted on various physical and o bservational 
groun ds to construct an a priori parameterisation (e.g. Brid le et ail 
l2003h . while others have attempted to reconstruct an u n- 
known spectrum directly from the data dWandl 19941: iB ridle et 
2003; Tocchini-Valentini, Douspis & Silk 2005; Hannestad 200s 
Shafieloo & Souradeed 120031) . In this paper we will mostly em- 
ploy the latter method, following closely the Bridle et al. analy- 
sis despite the obvious disadvantage of weakening constraints on 
remaining cosmological parameters and the creation of further de- 
generacies by the inclusion of more primordial parameters. One 
may, however employ a fully Bayesian extension of the Markov- 
Chain Monte Carlo (MCMC) method (see Sec. 12.1b to select the 
most appropriate model according to its Bayesian evidence (see 
ISlosar et alj|2003l: iBeltran et alj|2005[: iTrottal 120051 and most re- 
cently Mu kheriee. Parkinson & L iddle 2005). This method favours 
simple models with fewer parameters over a more elaborate model, 
unless the latter can provide a substantially better fit to the data. 

Using a parameterised spectrum requires constraining any 
analysis to a small sample of the huge number of possibilities, some 
of which have been outlined above. Observational features can re- 
fine the choice, for example the apparently low power observed in 
the first three CMB multipoles by WMAP. Although some authors 
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jEfstathioj 2003a) have suggested that these results are not statis- 
tically significant, in our analysis these data will always prefer a 
model with low power on these scales. We also test a broken spec- 
trum caused perhaps by a double field inflatio nary potential or mo - 
mentary pauses in the slow roll of the inflaton jBarriga et al.l200lh . 
lLasenbv & Doranl d2005h arrived at a spectrum naturally incorpo- 
rating an exponential fall off in power on large scales by consid- 
ering evolution of closed Universes out of a big bang singularity 
and with a novel boundary condition that restricts the total confor- 
mal time available in the Universe. The fact that this model fea- 
tures the type of cutoff that has been s uggested on phenomenologi- 
cal grounds from the data tefstathiou|[2 003b) make it an intriguing 
possibility. 



2 MODEL SELECTION FRAMEWORK 

2.1 Markov Chain Monte Carlo sampling 

A Bayesian analysis provides a coherent approach to estimating 
the values of the parameters, ©, and their errors and a method for 
determining which model, M, best describes the data, D. Bayes 
theorem states that 



P(O|D,A0 = 



P(D\&,M)P(&\M) 
P(D\M) ' 



(1) 



where P(0|D,M) is the posterior, P(D\&,M) the likelihood, 
P(&\M) the prior, P(D\M) the Bayesian evidence. Convention- 
ally, the result of a Bayesian parameter estimation is the posterior 
probability distribution given by the product of the likelihood and 
prior. In addition however, the posterior distribution may be used to 
evaluate the Bayesian evidence for the model under consideration. 

We will employ a MCMC sampling procedure to explore 
the posterior dis tribution using an ad apted version of the cos- 
moMC package i Lewis & Bridle 2002h with four CMB datasets ; 
WMAP dVerde et alJ2003l ; lHinshaw et alteOOSl i lkogut et alj2003h. 
ACBA R Ikuo et alj|2004h the latest V SA results jDickinson et all 
120041) and CBI dReadhead et alj|2004[). We a lso include the 2dF 
Galaxy Reds hift Surv ey (Perci val et al J 1200 If) , the Sloan Digital 
Sky Survey 1 Abazaii an et al. 2003) and the H ubble Space Tele- 
scope (HST) key project Freedman et al.l j200 lh . In addition to the 
primordial spectrum parameters, we parameterise each model using 
the following five cosmological parameters; the physical baryon 
density Qbh 2 ; the physical cold dark matter density Q. c h 2 ; the total 
energy density Qo (parameterised as Qk = 1 — fio); the Hubble pa- 
rameter h (Hq = h X lOOkms -1 ) and the redshift of re-ionisation 



2.2 Bayesian evidence and simulated annealing 

The Bayesian evidence can be defined as the average likelihood 
over the entire parameter space of the model: 



P(T>\&,M)P(&)d N &. 



(2) 



In general models with fewer parameters and hence small param- 
eter spaces will have larger evidences. Naively one could evaluate 
the evidence by sampling from the posterior distribution function 
randomly covering the entire parameter space and simply find the 
average likelihood value. In practise however the value of the like- 
lihood at its maximum is typically many orders of magnitude larger 
than that at any extremes of the space and so very many samples 



would need to be taken away from the peak in order for their in- 
clusion to make any appreciable difference. To overcome this one 
can use the numerical method of simulated annealing, also known 
as thermodynamic integration, to slowly increase the height of the 
peak in the likelihood relative to the surrounding background value. 
Using a random sampling of the space, such as MCMC, one can 
ensure that all of the parameter space, not just the area around the 
peak likelihood is explored. 

Using such an annealing schedule we draw samples from 
P(D|0, M) P(&\M) where A is the inverse temperature and is 
raised from w to 1. One can then define the evidence as a function 
of A as 



E(X) 



P{D\&,M) x P(&\M)d N &. 



(3) 



We require that the priors be normalised to unity over the parameter 
space considered thus E(X = 0) = 1 and we will also find it com- 
putationally simpler to work in the natural logarithm of evidences, 
therefore: 

InE(X) = lnP(0) + f ^^ dX 



dX 



1 dlnE(X) 



dX 



dX 



1 1 dE(X) 
E dX 



dX. 



Performing the derivative of E(X) gives 



ln£(A) 



1 / P(D\&,M) x \nP(D\&,M)P(&\M)d N & 
J P(Dj0, MY P(®\M)d N 



dX.(4) 



This gives the expectation value of In P(D [ 0, M) over the param- 
eter space defined by the extent of the priors. Thus we can simply 
sum over these values at each step AA in the annealing schedule 
and divide by the total number of samples to find the log evidence 
value 



ln.E(A) 



1 



E 

i=l 



lnP(D|0, M)i 



(5) 



The evaluation of this integral w as performed using a reverse an- 
nealing schedule (as set out in iBeltran et alj j2005h ) where A is 
lowered from 1 to once the Markov chain has found its station- 
ary point through a period of normal burn-in, which typically takes 
< 500 samples. The number of annealing steps N is not a chosen 
constant but is instead determined by the stopping criterion that the 
final log evidence value would not change by more than a given 
threshold (empirically set to 0.001). Accordingly A at each step is 



given as: 



X{N) = (1 - 6)' 



(6) 



where e is a user defined parameter of order 5 x 10 -5 . With this 
method log evidence values varying to within only one unit were 
successfully obtained from independent chains. The total number 
of sampler calls made during the evidence burn in was ~ 3000 de- 
liberately large to ensure the chain was at a stationary point. During 
the reverse annealing schedule typically ~ 15000 calls were made. 
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3 PRIMORDIAL POWER SPECTRUM 
PARAMETERISATION 

The correlation function £ over a range r of a density field, such as 
the matter density field in the Universe is defined as the product of 
the density contrast at two separate points, x and x + r 



£(r) = (5(x)5(x + r)). 



(7) 



The power spectrum P(k) is simply the inverse Fourier transform 
of the correlation function or the ensemble average power: 



P{k) = (\S k \ 2 ). 



(8) 



The basic inflationary paradigm inflates quantum zero-point oscil- 
lations of the inflaton to macroscopic scales, providing a source of 
density fluctuations which would then evolve via gravitational col- 
lapse to form the observed large-scale structure today. The power 
spectrum of these primordial fluctuations can then provide unique 
constraints on the dynamics of the inflationary epoch. The homo- 
geneity of the Universe on large scales suggests the simplest spec- 
trum and most obvious first candidate is scale invariant. This in- 
tuition is confirmed by the results of slow-roll inflation, predicting 
the slightly sloped power law spectrum 



P(k) = Ak" 



(9) 



where | n— 1 -C 1. It can be shown that the Fourier components 
are uncorrelated and have random phases meaning the power spec- 
trum encodes all of the stochastic properties of the density field. 



3.1 Harrison-Zel'dovich and Power-Law Parameterisations 

The simplest scale invariant power spectrum can be parameterised 
with one free parameter, an amplitude; A 3 as: 



P(k) = A s 



(10) 



This is known as the Harrison-Zel'dovich (H-Z) spectrum. Slow 
roll inflationary scenarios would be expected to imprint a slight 
slope, leading to the familiar power law form characterised by a 
spectral index n s (see Eqn.[9} A lin k between n s and the p otential 
driving inflation is well established dLvth and Riottoll 19981) but the 
spectrum is only an exact power law for an exponential inflationary 
potential, so in general the spectral index should be some function 
of scale n(k). Therefore we can also characterise the form of the 
spectrum as 



P(k) = A 



(feo) 



n-l + (l/2) ln(fe/fc )(dn/dln &+■■■) 



(ID 



where dn/d\n k is the running parameter n rll „. For most standard 
models we would expect n run ~ 0. The entire spectrum pivots 
about a scale ko which is also the position in fc-space at which the 
amplitude is set. In keeping with previous studies we have set this 
scale to O.OSMpc" 1 . 

The complete parameter estimates of our analysis using the 
data described above are shown in Table [T] The inclusion of the 
possibility of a non-flat geometry lowers the value of h consider- 
ably, a point to which we will return in Sec. 13.51 We find a best- 
fit single index spectrum (see Fig. [TJ that, despite the inclusion of 
more recent data and extension to non-flat geometri es, does not 
differ appreciably from the best-fit index found by Sperg el et al.l 
(2003). The inclusion of a running index (see Fig. [2} weakens the 
constraint slightly to 0.93 ± 0.05. To la level the constraint on 
n run provides only slight evidence for a dependence of n 3 on k. 
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Figure 1. Best fit single spectral index power law spectrum with lcr and 2<r 
errors (shown with shaded areas). [Note that, for clarity, the error in A s is 
not shown] 




k [Mpc 1 ] 



Figure 2. Best fit running spectral index power law spectrum with f a and 
2cr errors (shown with shaded areas). [Note that, for clarity, the errors in A s 
and n s are not shown] 



3.2 Large Scale Cutoff 

WMAP observations have pointed to lower than expected power 
on large scales. Although the statistical signific ance of these data 
points have been questioned (Efstathiou 2003b), currently models 
with a cutoff in power should be preferred. We did not attempt to 
model the form of the spectrum near the cutoff, instead parameter- 
ising the scale at which the power drops to zero, k c , with a prior of 
[0.0,0.0006] Mpc^ 1 : 



P(k) = 



0, 



k -5- kr 



(12) 



It is worth pointing out however that inherent cosmic variance lim- 
itations at this scale would make constraints of any form difficult to 
obtain. 

On small scales this spectrum behaves just as the single index 
power law and so constraints on the cosmological parameters re- 
main essentially unchanged. A cutoff is preferred (see Fig. [3j but 
the single index model is still not ruled out, illustrated by the non 
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Table 1. MCMC parameter constraints for H-Z, single index and running index parameterisations (mean ±1<t errors) 



Parameter 


H-Z 


Single Index 


Running Index 


Priors 


n s 

n b h 2 
n c h 2 

Q k 
li 

Zre 


27.1 ± 1.1 
0.0240 ± 0.0006 
0.116 ±0.009 
-0 026+ - 018 
63+ 05 
19.0 ± 2.6 


0.96 ± 0.02 

22.6 ± 1.1 
0.0229 ± 0.0009 
0.118 ±0.008 

-0 024+ ' 018 
u.uz^_ Qlg 

0.61 ± 0.05 
13.4 ±4.2 


0.93 ± 0.05 
-0.034 ± 0.032 

24.9 ± 1.2 
0.0225 ± 0.0012 
0.125 ±0.012 

-0 022+ 020 
u.uzz_ Q Qlg 

0.61 ± 0.05 
16.4 ±5.5 


[0.5,1.5] 
[-0.15,0.15] 
[14.9,54.6] 
[0.005,0.1] 
[0.01,0.99] 
[-0.25,0.25] 
[0.4,1] 
[4,30] 




0.0002 0.0004 

kJMpc 4 ] 



Figure 3. Marginalised histogram of the cutoff in power on large scales 



zero value of the likelihood function at k c — 0. The likelihood 
function peaks at a cutoff scale of 2.7 x 10~ 4 Mpc -1 and drops 
off markedly ther eafter toward s higher k c . This figure reproduces 
well the results of lBridle et al.l J2003h and despite a different math- 
e matical form for the c utoff gives a similar likelihood distribution 
to lNiarchou et alJd2004l) . it is also comfortab ly within the 2a upper 
limit of k c < 7.4 x 10" 4 Mpc -1 found bv lContaldi etail <2003h 
using an exponential cutoff . 

3.3 Broken Power Spectrum 

Rapid cooling of the Universe during inflation can result in spon- 
taneous symmetry breaking phase transitions which could interrupt 
the potential driving inflation leading to one or more sudden de - 
partures in scale invariance lasting ~ le fold tearriga et alj|200lh . 
When this occurs, the mass of the inflaton field changes suddenly, 
thus interrupting inflation. Although predicting the form of this 
interruption is not trivial, it is clear we should observe a drop in 
power, i.e. a break in the primordial power spectrum, after which 
the spectrum should return to scale invariance. We will assume the 
spectrum has the general form 

{A, k k s 

Ck a ~\ k a <k^k e (13) 
B, k > fc e 

where the values of C and a are chosen to ensure continuity. Four 
power spectrum parameters were varied in this model: the ratio of 
amplitudes before and after the break A/B with prior [0.3,7.2]; 
k s indicating the start of the break with prior [0.01, 0.1] Mpc -1 ; 
\n(k e /k s ) to constrain the length of the break with prior [0, 4] and 
normalisation A with prior [14.9, 54.6] x 10 -8 . We also placed a 
prior that k e could not exceed 0.1 Mpc -1 . 




0.02 0.04 0.06 1 1.2 1.4 




1% A/B ln(k e /k s ) 

Figure 4. Marginalised ID and 2D probability constraints on k s ln(fc e /fe 3 ) 
and A/B, 2D constraints plotted with 1<t and 2a confidence contours. 

Constraints on these parameters are shown in Fig. [4] in very 
good agreement with the Bridle et al. analysis showing that this 
parameterisation is robust in the extension to closed models. From 
the preference for a spectral index lower than unity in Sec 13.11 it 
is clear that a drop in power would be preferred by the data, es- 
pecially if the transition was smooth and extended as this would 
mimic a tilted spectrum. Indeed this effect is seen by the large peak 
in the likelihood surface at very low values of ln(fc e /fc s ) (i.e. a large 
difference in k e and k s implying an extended break). More inter- 
esting however is the preference for a sudden drop on large scales 
k ~ 0.025Mpc _1 , this could be evidence for a phase transition in 
the early Universe, or as Bridle et al. point out, could be an artifact 
of the combination of WMAP and 2dF datasets. Example spectra 
featuring both extended and sharp breaks are shown in Fig. [5] 



3.4 Power Spectrum Reconstruction 

The data that can now be used to constrain the primordial power 
spectrum are generally analysed in the framework of a specific 
model, as we have done up to now. It is possible however that cer- 
tain features of the observational data are being overlooked with 
this method. In order to reveal any such features we have divided 
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Figure 5. Examples of two broken spectra; with an extended and sudden 
transition in power. 
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Figure 6. Reconstruction of the primordial power spectrum in 8 bands of k, 
compared to the best fit H-Z spectrum (dotted line). 



the spectrum into eight bins in wavenumber k and allowed the am- 
plitude in each bin to vary over the same range as the single am- 
plitude As for the power law spectrum. Instead of using simple 
'topha t' bin amplitudes we used the same approach as lBridle et al.l 
d2003h . linearly interpolating between the bins in k. Since cur- 
rent interest has concentrated on large scale departures from scale 
invariance we confined the study to bins, of amplitude a n , be- 
tween 0.0001 and 0.11 Mpc -1 parameterised logarithmically with 
fcj+i = 2.75ki so that 



P(k) = 



ki < k < fci+i 
k ^ kn 



(14) 



We show in Fig. [6] the reconstructed spectrum relative to the 
best fit Harrison-Zel'dovich (H-Z). The effect of degeneracies be- 
tween neighbouring bins is minimised with the use of cosmoMC's 
new parameterisation to explore the A s — r degeneracy, however 
there is still considerable variability at the la level. Our results 
are mostly in agreement with a H-Z spectrum, altho ugh we also 
find th e feature observed at around 0.06 Mpc -1 by iBridle et al.l 
d2003l) . We observe a slight decrement in t he power on scale s be- 
low 0.001 Mpc -1 , confirmi ng results by Hannestad) d2004l) and 
IShafieloo & Souradeed J2003h . this effect would produce a similar 
decrement in the CMB spectrum at low I. 



3.5 Closed Universe Inflation 

lLasenbv & Doranl d2005l) arrived at a novel model spectrum by con- 
sidering a boundary condition that restricts the total conformal time 
available in the Universe, and requires a closed geometry. The re- 
sultant predicted perturbation spectrum encouragingly contains an 
exponential cutof f (as previously suggested phenomenologically by 
lEfstathioul l2003a) at low k which yields a corresponding deficit in 
power in the CMB power spectrum. The shape of the derived spec- 
trum was parameterised by the function: 



P(k) = A(l - 0.023y) 2 (l - exp(-(y + 0.93)/0.47)) 2 



x 3 x W) > -0.93 



(15) 



where J/ = ln (jWToo 

Unfortunately it has not been possible to produce a routine 
that can recalculate the form of this spectrum (see Fig. [7} in non- 
primordial parameter space fast enough to perform a full MCMC 
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Figure 7. Lasenby & Doran spectrum (short-dashed) shown with best fitting 
H-Z (dot-dashed), single index with a cutoff (dotted), without (full) and a 
running index (long-dashed). 



analysis. The authors did, however suggest a best-fit cosmology 
(D. b h 2 = 0.0224, Q, cdm h 2 = 0.110, h = 0.6 and a total energy 
density fin = 1.04) from which we could fit the primordial spec- 
trum using just an amplitude A with prior [14.9, 54.6] x 10 -8 . 
Although such a low Hu bble parameter wou ld have been ruled out 
from previous analyses dSpergel et ai]|2003l) the extension to non- 
flat geometries weakens the constraint considerably, as discussed 
earlier. The best fitting spectrum (A = 29.82 ± 0.19) is shown with 
the H-Z, single index, cutoff and running index models for compar- 
ison in Fig|7] Encouragingly the 'turn-over' scale of the L+D spec- 
trum is close to that found in the reconstructed spectrum at about 
0.001 Mpc -1 , while at large k the spectrum successfully mimics a 
power law spectrum with spectral index n s » 0.96. 



4 APPLICATION TO SIMULATED DATA 

In order to test the model selection algorithm we produced a se- 
ries of simulated data sets with known cosmological parameters 
and a particular primordial power spectrum parameterisation. We 
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Figure 8. Simulated cosmic variance limited data to I = 2000 (convolved 
with Gaussian noise) and the model from which it was produced (dashed 
line). 



Table 2. Differences of log evidences for four parameterisations using cos- 
mic variance limited data. 



Model 


lnE A - InE 


Constant n 


0.0 ± 2.4 


Running 


-4.7 ± 2.2 


Cutoff 


-1.1 ± 2.9 


Broken 


-36.7 ±3.3 



produced two sets: the first intended to represent the quality ex- 
pected from the forthcoming Planck satellite mission, the sec- 
ond a combination of simulated WMAP first year CMB data and 
Sloan Digital Sky Survey (SDSS) LSS data. At least eight separate 
Markov chains were produced giving eight estimates of the evi- 
dence and an associated variance. The chosen model was flat with 
n b h 2 = 0.022, n cdm h 2 = 0.120, h = 0.72, optical depth to re 
ionisation of r = 0.15 and with primordial parameters n 3 = 0.97 
wdA 3 = 2.5 x 10" 9 . 



4.1 Simulated Planck Data 

We approximated mock Planck data as being cosmic variance lim- 
ited up to I of 2000, which was sufficiently accurate to determine 
the spectral index n s in such a parameterisation to an accuracy 
of better than 2% and was produced using C;s from the power 
spectrum generator CAMB (l iLewis. Cha llinor & Lasenby 200(3)) to 
which we added Gaussian noise with a cosmic variance standard 
deviation of 

2 



-C ■ 



(16) 



1(1 + 1) 

The CMB spectrum of the chosen model is shown in Fig. [8] along 
with the data points used including simulated noise. Parameter con- 
straints confirm that the sampler is capable of extracting this model 
from the dataset (see Fig. [9}. 

Since the absolute evidence value can only be compared be- 
tween models using the same dataset it is conventional to quote 
ratios with respect to a reference model, which in this case will be 
the single spectral index. Our results (shown in Table [2} success- 
fully select the correct model, although with high variance between 
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Figure 9. Recovered MCMC marginalised parameter constraints using cos- 
mic variance limited data. 



chains. In particular, the cutoff model is, within estimated la er- 
ror, a very close second. This is not too surprising given that the 
cutoff spectrum is simply a single index spectrum above the appro- 
priate cutoff wave vector fc c , at which scale cosmic variance is large 
enough to prevent one obtaining strong constraints. The question of 
what difference in log evidence one requires to say with some cer- 
tainty which model the data prefers is not a trivial one and depends 
to some exten d on 'intuition'. A useful guide has been given by 
Ijeffrevsl Jl96lh where a log evidence difference: AlnE < 1 is not 
significant, 1 < AlnE < 2.5 significant, 2.5 < Aln_E < 5 strong 
and AlnE > 5 decisive. Using this criterion we can correctly rule 
out, with some confidence, the running index spectrum and com- 
pletely discount the broken model. The combination of some LSS 
data sets would doubtless allow a distinction between the cutoff 
and single index spectrum though this test has not, as yet been per- 
formed. 



4.2 Simulated WMAP & SDSS Data 

Two further datasets were created to mimic the current data used in 
our full analysis: a simulated first year WMAP CMB power spec- 
trum and large-scale structure (LSS) data similar to the SDSS mat- 
ter power spectrum. Both spectra were again produced using the 
CAMB generator using the same single index model chosen above. 
The mock WMAP Cis were created by adding Gaussian noise us- 
ing the calculated WMAP errors at multipoles between 2 and 856 
corresponding to the binning used by the WMAP team (see Fig. 
UOt . SDSS data was produced (see Fig.ll It using the binning in k 
space and calculated errors used by (?). As is well known the com- 
bination of LSS and CMB data is capable of breaking a number of 
degeneracies in parameter space and can produce tighter constraints 
than with CMB data alone (see Fig. 1121 ). Although differences in 
the evidence estimates made with this realistic data are somewhat 
smaller than with simulated Planck data it is still possible to deter- 
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Figure 11. Simulated SDSS data and the matter power spectrum from 
which it was produced (dashed line) 



mine the correct model, though not decisively (Table [5} according 
to Jeffreys criteria. 



5 APPLICATION TO REAL DATA 

Given the success of model recovery using simulated data, we can 
now turn back to real data with some confidence. All of the param- 
eter estimations in Sec.[3]were accompanied by at least eight sepa- 
rate log evidence estimates. We can divide this part of the analysis 
into two main sections: firstly a full parameter space exploration for 



Table 3. Differences of log evidences of four models using a combination 
of simulated WMAP and SDSS data. 



Model 


ln.E A - ln_E 


Constant n 


0.0 ±0.8 


Running 


-2.3 ± 0.8 


Cutoff 


-1.2 ±0.6 


Broken 


-1.1 ± 1.5 



3 3 4 

i°g[io 10 Aj 



Figure 12. Recovered MCMC marginalised parameter constraints for the 
chosen model using simulated first year WMAP + SDSS data. 



the H-Z, single index, running, cutoff, broken and binned models; 
and secondly estimations made only with the reduced primordial 
parameter space (i.e. only those parameters affecting the form of 
the initial power spectrum) within a predetermined cosmology, to 
include the Lasenby & Doran model. This division was necessary 
because of the restrictions in the form of the Lasenby & Doran 
spectrum discussed in Sec. 13.51 



5.1 Full cosmological parameter space exploration 

Statistical uncertainty dominates the results of the full parameter 
space study (see Table |4), leaving it difficult to draw any mean- 
ingful conclusions. All of the central mean values, with the excep- 
tion of the binned spectrum, are found to differ by less than a unit 
in log evidence, which, according to Jeffreys criteria, provides an 
insignificant difference in evidences. However ranking these mean 
values we find that the cutoff spectrum is marginally pr eferred over- 
all. Al though a better chi-squared fit was found by ISpergel et al.l 
(2003) for a running parameterisation we find no conclusive evi- 
dence in its favour in this study. Indeed it is worth noting, any pos- 
sible preference for a typical running spectrum(?i rtln ~ —0.035) 
could possibly be due to its large scale power decrement rather 
than any running. The reconstructed spectrum, revealing any possi- 
ble structure would be expected to provide the best possible model 
fit -however the large increase in parameter space has clearly had 
a detrimental effect producing a 'significant' evidence result dis- 
favouring the model. In order to reduce the statistical uncertainty 
in these estimates to the point of distinguishing between the best 
fit models we would require of order 50 seperate runs for each 
parameterisation which at present is not computationally feasible. 
However the use of the alternative method of nested sampling (see 
Mukh eriee et ai] |2005h could allow the reduction of uncertainties 
to the required level in a future study. 
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Table 4. Differences of log evidences with respect to single index model 
using the full cosmological parameter space. 



Model 


\nE h - InE 


Constant n 


0.0 ± 0.7 


H-Z 


0.4 ± 0.8 


Running 


0.5 ± 0.7 


Cutoff 


0.8 ±0.8 


Broken 


-0.3 ± 0.8 


Binned 


-1.8 ±0.7 



Table 5. Differences of log evidences (for primordial parameters) for all 
models with respect to single index model within cosmology: f2o = 
1.04, Q b h 2 = 0.0224, H = 60, fl cdm h 2 = 0.110. 



Model 


ln£ A - ln-E 


Constant n 


0.0 ± 0.6 


H-Z 


-0.4 ± 0.5 


Running 


-2.1 ±0.5 


Cutoff 


0.2 ± 0.6 


Broken 


-0.5 ± 0.7 


Binned 


-5.2 ±0.7 


Lasenby & Doran 


0.9 ± 0.6 



5.2 Primordial parameter space exploration 

The Lasenby & Doran spectrum had previously been calculated 
for a cosmological model described by Qn = 1-04, Qbh 2 — 
0.0224, h = 0.6, Qcdmh 2 = 0.110, and provided a good fit to 
WMAP and higher resolution data. We therefore used this cos- 
mology as a starting point. Again the results (see Table [5} have 
been normalised to the single index spectrum, and the similarity 
in evidences make it difficult to distinguish between the H-Z, sin- 
gle index and cutoff. We can however clearly rule out the running 
spectrum in this case with some certainty. Though not a convinc- 
ing result in this case, the mean evidence of the Lasenby & Doran 
spectrum lies slightly higher , quan titatively confirming the 'hand' 
fit used in lLasenbv & Doranl J2005t) . 

With this relative success for a restricted cosmology, it made 
sense to consider the evidence rankings within the best fit con- 
cordance cosmology, in this case the parameter estimates found 
from the full MCMC simulations of the power spectrum. A general 
package for predicting the power spectrum for arbitrary parameters 
within the Lasenby & Doran model is not currently available, how- 
ever we know the form to be fairly stable with varying cosmological 
parameters. Hence we used the Qq = 1.04 model described above 
as a template to compare with the results from the full MCMC fit- 
ting, which gave parameters £1q = 1.024, flbh 2 = 0.0229, h = 
0.61, fl c dmh 2 = 0.118. Although these correspond to a different 
cosmology from that used for the Lasenby & Doran case, the cos- 
mological parameters do not deviate greatly from these across the 
different parameterisions, in particular the cutoff and running cases, 
hence we feel an evidence comparison is justified. This result (Ta- 
ble^, shows the L+D model to lie roughly 3 units of log evidence 
ahead of the cutoff case, yielding an evidence ratio of roughly 20, 
which according to the Jeffreys ranking system provides a 'strong' 
model selection. 



Table 6. Differences of log evidences (for primordial parameters) for all 
models with respect to single index model within the current concordance 
cosmology: U = 1.024, Q b h 2 = 0.0229, h = 0.61, VL cdm h 2 = 0.118, 
as compared to the Lasenby & Doran model (with parameters as in the 
previous table) 



Model 


IilEa - to.E 


Constant n 


0.0 ± 0.5 


H-Z 


-4.4 ± 0.5 


Running 


-0.8 ± 0.6 


Cutoff 


0.4 ± 0.5 


Broken 


-2.7 ± 0.6 


Binned 


-6.1 ±0.6 


Lasenby & Doran 


4.1 ±0.5 



6 CONCLUSIONS 

We have performed a full Bayesian analysis of various parameter- 
isations of the primordial power spectrum which includes not only 
the estimation of cosmological and spectral parameters but also a 
value of the Bayesian evidence. This method of model selection 
was tested using simulated current and future data, though it is dif- 
ficult to make conclusive determinations in most cases using real 
data at present. Consistently high evidence values were obtained 
for those models incorporating a reduction in power at low k, a nat- 
ural result given the 1 st year WMAP data. In particular we found 
strong evidence in favour of the Lasenby & Doran spectrum within 
a limited MCMC analysis. 
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